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ABSTRACT 

The physical nature of thermal composite supernova remnants (SNRs) remains controversial. 
We have revisited the archival XMM-Newton and Chandra data of the thermal composite SNR 
Kesteven 41 (Kes 41 or G337.8—0.1) and performed a millimeter observation towards this source 
in ^^CO, ^^CO, and C^®0 lines. The X-ray emission, mainly concentrated toward the southwest¬ 
ern part of the SNR, is characterized by distinct S and Ar He-like lines in the spectra. The X-ray 
spectra can be fitted with an absorbed non-equilibrium ionization collisional plasma model at a 
temperature 1.3-2.6 keV and an ionization timescale 0.1-1.2 x 10^^cm“^s. The metal species 
S and Ar are over-abundant, with 1.2-2.7 and 1.3-3.8 solar abundances, respectively, which 
strongly indicate the presence of a substantial ejecta component in the X-ray emitting plasma of 
this SNR. Kes 41 is found to be associated with a giant molecular cloud (MC) at a systemic LSR 
velocity —50kms“^ and confined in a cavity, delineated by a northern molecular shell, a western 
concave MC that features a discernible shell, and an H I cloud seen toward the southeast of the 
SNR. The birth of the SNR in a pre-existing molecular cavity implies a mass > I8M0 for the 
progenitor if it was not in a binary system. Thermal conduction and cloudlet evaporation seem to 
be feasible mechanisms to interpret the X-ray thermal composite morphology, while the scenario 
of gas-reheating by the shock reflected from the cavity wall is quantitatively consistent with the 
observations. An updated list of thermal composite SNRs is also presented in this paper. 

Subject headings: ISM: individual (G337.8—0.1 = Kesteven 41), — ISM: molecules, — supernova rem¬ 
nants 


1. Introduction 

Massive stars rapidly evolve to the end of their 
nuclear-burning lifetimes and explode as core¬ 
collapse (CC) supernovae before they can move 
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far away from where they formed. Dozens of CC 
supernova remnants (SNRs) have been found to 
be interacting wi th molecular clouds (MCs) (see 
Jiang et al. 1 1201(1 and references therein). Such 


SNRs are a crucial probe for hadronic interactions 
between the shock accelerated protons and the as¬ 
sociated MCs, and many detections of them have 
been made based on 7-ray observations in the cur¬ 
rent Fermi and HESS era. Among the m are a num¬ 
ber of so-called th e rmal composite ( Jones et al.l 
19981: IWilner et al.l Il998^ or mixed-morphology 


remnants (|Rho fc Petrelll998l) : this class of SNRs 
is characterized by the coexistence of radio shells 
and centrally brightened thermal X-ray emission. 
The large majority of thermal composite SNRs 
have been found to be interacting with adjacen t 


MCs ( Creen et al.lll997HYusef-Zadeh et al.ll2003 ). 
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Unlike the regular composite type SNRs in which 
the central X-ray emission is clearly known to be 
powered by pulsars, the nature of the X-ray kernel 
in the thermal composites remains a controversial 
issue. 


Several attempts have been made to interpret 
the X-ray morphology of the thermal composites 
and some mechanisms to explain the center-filled 
X-ray morphology have been proposed, including 
hot interior with radiatively cooled rim, thermal 
conduction in the interior hot gas, gas evapora¬ 
tion from the shock-engulfed cloudlets, projection 
effects of the one-sided shock-cloud interaction, 
metal enrichment in the interior, shock reflected 
inward from wind-blown cavity wall, amonst 
other proposed m echanisms (see a summary in 
Chen et al.l l2008l and references therein). More 
studies of thermal composite SNRs are needed to 
test these mechanisms. For this purpose, we here 
investigate the properties of the interior hot gas 
and the environs of SNR Kesteven 41 (Kes 41 for 
the remainder of this paper). 

Kes 41 (G337.8—0.1), a southern-sky SNR first 
discovered with the Molonglo Obse rvatory Synthe- 
sis telescope (MOST) at 408 MHz ( Shaver fc Gosd 
1970ll . is shown to be centrally brightened in X- 
rays within a distorted radio shell by an XMM- 
Newton observation and i s therefore classifie d as 
a thermal composite SNR ( Combi et al.l 2008 1. In 
a spectral analysis, the X-ray emission was sug¬ 
gested to arise from a hot (~ 1.4 keV) gas of nor¬ 
mal metal abundance; however, only the MOS 
dat a of the XMM-New ton have been analyzed so 
far (jCombi et al.ll2008h . Presently, the available 
X-ray observations allow a more thorough analy¬ 
sis of the physical properties of Kes 41. 

Moreover, Kes 41 has been found to be in¬ 
teracting with an adjacent MC as indicated by 
the 1720 MHz hydroxyl radical ( OH) maser emis¬ 
sion detected on th e radio shell ( Koraleskv et al.l 
19981 : Caswelll[2004ll . The OH satellite line masers 
at 1720 MHz are widely a ccepted as a sign¬ 
post of SNR-MC interaction (Lockett et al. 1999t 


Frail fc Mitchelll 1998HWardle fc Yusef-Zadehl200^ . 

The SNR-MC interaction may probably con¬ 
tribute to the high-en ergy y-ray excess near Kes 41 
detected by EGRET (ICasandiian fc Grenieijl2008 1 
and Eermi i Ergin fc Ercan 2012ll satellites. The 
distance to the remnant is estimated to be 
12.3 kpc from the OH maser detection at the lo- 


cal standard of rest ( LSR) velocity —45 kms“^ 
( Koraleskv et al. 1998h and the H 1 absorption 
that places it beyon d the tangent point at 7.9 kpc 
(ICaswell et al.lll975ll . As a follow-up study of the 
association of Kes 41 with an adjacent MC, we 
have performed a millimeter observation toward 
Kes 41 of CO transition lines. 


In this paper we carry out an X-ray analysis 
of SNR Kes 41 with archival XMM-Newton (MOS 
and PN) and Chandra observation data and inves¬ 
tigate the ambient interstellar environment with 
our new CO line observation as well as an archival 
H I survey. 


2. Observations and Data Reduction 


2.1. XMM-Newton and Chandra X-ray Ob¬ 
servations 


The XMM-Newton X-ray satellite observed 
Kes 41 on September 26th 2005 using both 


the EPIC-MOS and the EPIC-PN (Combi et al 


, centered at (RA (J2000.0) = 16^39“06!5, 
DEC (J2000.0) = -46°57'58"0). The raw EPIC 


data are cleaned with the latest version of the 
Standard Analysis System (SAS-V13. 0 . 3 ) 0 , of 
which the Extended Source Analysis Software 
(ESAS) package is mainly used because Kes 41 
appears as an extended source in these observa¬ 
tions. After removing time intervals with heavy 
proton darings, the edective net exposure time for 
both the EPIC-MOS and the EPIC-PN is 35 ksec. 
We exclude dux from all the detected point sources 
before performing spatial and spectral analyses. 


Chandra observed Kes 41 with the Advanced 
CCD Imaging Spectrometer (ACIS) on June 2011 
incidentally in an observation towards high-mass 
X-ray binaries in the Norma region (PI: John Tom- 
sick). There are six sections of the observation 
that cover Kes 41 region (see Table [T] for observa¬ 
tional parameters). We reprocess the events dies 
using the Chandra Interactive Analysis of Obser¬ 
vations (CIAO) data processing software (Version 
4 . 5 ) 0 . In the reprocessing, bad grades are dltered 
out and good time intervals are reserved. Flux 
from all the detected point sources in the held of 
view (FOV) are also removed. 


^See http://xmm.esac.esa.int/sas/ 
^See http://cxc.harvard.edu/ciao/ 
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2.2. CO Observation and Data 


3. Analyses and Resnlts 


The spectroscopic mapping observation of CO 
molecular lines, i.e. (J=l-0) at 115.271 GHz, 
13CO (J=l-0) at 110.201GHz, and C^^O (J=l- 
0) at 109.782 GHz, was made with the Mopra 
22-m telescope in Australia on 7-11 July 2011. 
We did on-the-fly (OTF) mapping with the field 
of ^ 11' X 10' covering the entire angular ex¬ 
tent of Kes 41. The 3-mm band receiver and the 
Mopra spectrometer backend system were used in 
zoom mode. The bandwidth in this mode is about 
137.5 MHz (corresponding to about 360kms“^ for 
^^CO (J=l-0)) with a spectral resolution of about 
33.6 kHz (corresponding to an LSR velocity reso¬ 
lution of ^ 0.09kms“^). During the observation 
the system temperature was around 280 K and the 
pointing accuracy was mostly better than ~ 3.5". 
The half-power beam size is ~ 30" at 115 GHz, 
with the ma in beam efficiency ~ 0.42 at the same 
wave band ( Ladd et’afil 200511 . In this paper, all 


the images are presented in antenna temperature, 
but radiative temperature after calibration is used 
to derive the physical parameters of the molecular 
gas. 

The software Livedata/Gridzill^I is used 
to calibrate the OTF-mapping data and combine 
them into three-dimensional data cubes, with two 
spatial dimensions (longitude/latitude) and one 
spectral dimension (LSR velocity). The subse¬ 
quent data reduction, such as subtracting fitted 
polynomial baselines, are performed using our own 
script written in Python. Three data cubes of 
diverse GO lines are yielded with grid spacing of 
30" and velocity resolution of ~ 0.1 kms“^. Fi¬ 
nally, the mean rms noise levels of the main beam 
brightness temperature are about 0.3, 0.15, and 
0.15 K for the ^^GO (J=l-0), ^^GO (J=l-0), and 
G^®0 (J=1~0) lines respectively. 


2.3. Other data 

The archival radio data from the MOST Galac¬ 
tic plane survey at 843 MHz (with an angula r res¬ 
olution of ^ 43"; Whiteoak fc Greei] Il99fi l. and 
the SGPS H I data (with an angular resolution 
of ^ 2' and a spectral resolu tion of ~ 0.8kms“^; 
McGlure-Griffiths et al.ll200^ are also used in our 
multi-wavelength study of Kes 41. 


3.1. X-ray properties 

3.1.1. Spatial analysis 

In Figure [T] we present a broadband (2.0- 
7.2 keV) X-ray image of Kes 41 made from the 
EPIG-MOS and EPIG-PN data, as coded in blue. 
In this image, we have removed the bad pixels, 
masked the detected point sources for both MOS 
and PN, merged the MOS and PN imaging data 
with the quiescent particle background subtracted 
and exposure map correction applied, and adap¬ 
tively smoothed the image (using the tool fadapt of 
HEASOFtH, version 6.14). It can be clearly seen 
that the X-ray emission is internally filled with¬ 
out limb brightening as contrasted to the shell-like 
feature in radio (shown in white contours). The 
centroid of the X-ray emission sits in the southwest 
(SW) region of the northeast-southwest elongated 
extent of the SNR. 

In Figure 2 we present X-ray images that are 
produced in the 2.0-2.8 keV, 2.8-3.5 keV and 3.5- 
7.2 keV energy ranges, respectively, following the 
similar procedure for the broad band image. These 
images illustrate the photon energy dependence 
for broadband intensities. The narrow bands cor¬ 
respond to the sulfur line, argon line (see S I3.1.2p . 
and higher energy emission, respectively. In both 
the soft (2.0-2.8 keV) and hard (3.5-7.2 keV) im¬ 
ages, there are two X-ray brightness peaks located 
in the geometric center and the southwestern re¬ 
gion. The hard X-ray image looks relatively bright 
in the SW compared with that in the NE. 

All six short sections of the Chandra observa¬ 
tion described in Table [T] are also used to gener¬ 
ate a merged broadband image of Kes 41, but the 
location of SNR is near to the edge of the GGD 
chips in each section of observation, resulting in a 
decrease of the angular resolution. The Chandra 
X-ray image of the SNR shows no significant dif¬ 
ference from the XMM-Newton image and is not 
presented here. 

3.1.2. Spectral analysis 

Figure [3] shows the simultaneous fit of the 
XMM-Newton and Chandra spectra of Kes 41. 
The spectral extraction regions for the diffuse 


^See http://www.atnf.csiro.au/computing/software/livedata/ 


^See http://heasarc.gsfc.nasa.gov/lheasoft/ 


3 














emission source and the local background are de¬ 
fined in Figure |4] with all the detected point 
sources removed. The spectral extraction is such 
as to include the brightest portion of the diffuse X- 
ray emission, which is located in the southwestern 
part of the SNR and is within the FOVs of all the 
XMM-Newton MOS and PN and Chandra-AClS 
detectors. Because of the insufficiency of event 
counts in each spectrum of the XMM-Newton 
EPIC-MOS and Chandra observation, we merge 
the MOSl and MOS2 spectra using the HEA- 
SOFT tool addascaspec, and merge three Chan¬ 
dra spectra (see Table [I]) using combine spectra in 
CIAO 4.5 to increase the statistical quality. The 
other three sections of the Chandra observation 
hardly cover the main area of the X-ray emis¬ 
sion and are thus not considered. The spectra are 
all adaptively regrouped to achieve a background- 
subtracted signal-to-noise ratio of 3 per bin. In 
these spectra, there are apparent emission line fea¬ 
tures of metal elements S (~ 2.43 keV) and Ar 
(~ 3.16keV), as well as a hint of Ca (~ 3.9keV), 
which confirms the thermal origin of the X-ray 
emission. There is no emission below ~ 2.0 keV, 
which indicates a large value of line of sight (LOS) 
hydrogen column density Nn. 


The XMM-Newton and Chandra spectra can 
be jointly fitted with an absorbed (using XSPEC 
model phabs) single temperature non-equilibrium 
ionization (NEI) model (vnei, version 2.0) and the 
results are summarized in Table [5] In the spec¬ 
tral fit, we allow the abundances of S and Ar, 
which have distinct line features, to vary and set 
the abun dances of the other metal elements equal 
to solar ( Anders fc Greyest 1989l l. The thermal 
X-ray emitting gas is at a temperature ~ 1.3- 
2.6 keV with an ionization parameter Uet ^ 0.1- 
1.2 X 10^^cm“^s. The best fit abundances of the 
two elements are moderately elevated (^ 1.2-2.7 
and ^ 1.3-3.8 times solar, respectively), indicat¬ 
ing that the hot gas is metal-enriched. 

To see where the metal-rich materials’ emission 
arises, we present the equivalent width (EW) im¬ 
ages of S and Ar lines in Figure [5l These images 
are con structed using a met hod similar to that 
used in Jiang fc Chenl ( 2010( 1. so as to maximize 
the statistical quality. Assuming that the line- 
to-continuum ratio is intrinsically invariable for 
various detectors, we sum the X-ray counts from 
XMM-Newton EPIC-MOS and PN and Chan¬ 


dra ACIS observations band-by-band according 
to componendo. We then rebin the data using 
an adaptive mesh, with each bin in each narrow 
band image containing about 10 counts. The in¬ 
tensity of the line emission with the interpolated 
continuum component subtracted is used as the 
numerator of the ratio, and the background (es¬ 
timated from the region surrounding the SNR) is 
subtracted from the continuum intensity which is 
used as the denominator. Figure [S] shows that the 
EW values of both S and Ar lines are only signifi¬ 
cant in the interior of the southwestern part of the 
SNR. The caveat here is that the derived contin¬ 
uum images might have lower statistical qualities 
than the corresponding line emission images and 
not all the bins can reach 3a significance especially 
outside the X-ray emission region. 

3.2. The Interstellar Environment 

3.2.1. Spatial distribution of the ambient clouds 

The average CO spectra over the FOV are 
shown in Figure |6l There are a number of spectral 
components along the LOS at the LSR velocities 
from WsR = —130 kms“^ to Ksr = 0 kms“^. 
The LSR velocity of the 1720 MHz maser spot, 
^ —45 kms“^, is in the right wing of the line 
component peaked at ~ —50 kms“^ which is com¬ 
monly revealed in the ^^CO (J=l-0), ^^CO (J=l- 
0), and C^®0 (J=l-0) spectra (see Figure (S]). 
By an examination of the channel maps of the 
~ —50kms“^ line component, we find a morpho¬ 
logical correspondence between the western radio 
boundary of the SNR and a concave surface of 
MC seen in the ^^CO emission at velocity inter¬ 
val ^ —55—52 kms“^ (see Figure [T]). The con¬ 
caved MC is also remarkable in the integrated 
^^CO intensity map in interval —70—40kms“^ 
within the velocity range of the line profile (see 
Figure [1}. In ^^CO emission, a bow-like structure 
at ~ —50—48kms“^ seems to match the west¬ 
ern radio shell; furthermore, a molecular shell at 
~ —61—58kms“^ perfectly follows the western 
radio shell of the SNR (see Figure (H). In addi¬ 
tion, a section of molecular shell in ^^CO emis¬ 
sion at —47.5kms“^ is seen to closely follow the 
northern radio shell (see Figure 0. Such corre¬ 
spondence is actually only seen in this spectral 
component by examining the data cubes across 
the entire LSR velocity span. Around the OH 
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maser position, we can also see some dense molec¬ 
ular gas in the intensity maps (Figure |9|) of C^®0 
emission, which is optically thin and traces molec¬ 
ular cores. The combination of the OH maser and 
the spatial features of CO emission presented here 
clearly demonstrates that SNR Kes 41 is in phys¬ 
ical contact with the MC at a systemic velocity 
Vlsr ^ —50kms“^ in the west and the north, 
although the broadened CO line profiles of the 
~ — 50kms“^ molecular component due to shock 
disturbance can not be discerned because of line 
crowding along the LOS. 

It is noteworthy that a dense H I cloud at 
Vlsr ~ — 52kms“^ is situated to the southeast 
of the SNR (see Figure [I] where the H I inten¬ 
sity image is presented in velocity interval —55- 
—50 kms“^), and the southeastern radio shell of 
the SNR seems to follow the surface of the H I 
cloud. SNR Kes 41 appears to be confined in a 
cavity enclosed by the molecular cloud in the west 
and northwest and the H I cloud in the southeast. 
The archival data show that this H I cloud has an 
angular size of ~ 18'. 


the second and the third method, we have as¬ 
sumed a local thermodynamic equilibrium and 
^^CO (J=l-0) line being optically thick, with an 
excitation temperature of 12.9 K (as obtained from 
the line peak of (J=l- 0 )). 

Prior to deriving the column densities of 
and C^® 0 , we have derived the optical depths 
t(^^CO) « 0.6 and t(C^®0) « 0.1, and therefore 
we treat the and C ^®0 emission as optically 

thin. The A^(H 2 ) value obtained from the third 
method, ^ 5.1 x 10^^ cm“^, is a little greater than 
those obtained from the other two methods, and 
it is probably because C ^®0 traces denser clouds. 


The archival data show that the H I cloud 
southeast to Kes 41 also has a spectral peak at 
—50kms“^ and similar velocity span to the spec¬ 
tral profile of the —50kms“^ component. 

The column density of H I gas (Tg ~ 110 K) is 
estimated to be ^ 5.1 x lO^^cm”^, using the re¬ 
lation iV(H I)=1.823x lO^^TgAti on the op tically 
thin assumption ( Dickev fc Lockmanlll99C 


4. Discussion 


3.2.2. Parameters of the ambient clouds 


4.1. SNR environment 


We estimate the gas parameters of the —50 km s“^ 
molecular component in velocity range -70- 
—40 kms“^ with which the SNR is interacting, 
and the results are listed in Table El The col¬ 
umn density and mass of the molecular gas de¬ 
rived may be an upper limit, because there may 
be irrelevant contributions in the complicated 
line profiles in and emission. For 

the estimate, we adopt a box region including 
the main part of the MC in the FOV (see the 
—53.5kms“^ panel of Figure 111). Three meth¬ 
ods have been utilized to estimate the average 
molecular column density A(H 2 ), which give a 
similar value iV(H 2 ) ~ 3-5 x 10^^cm“^. In 
the first method, we use the mean CO-to-H 2 
mass conversion factor (kn own as the X-fac tor) 
1.8x 10^° cm“^ K“^ km“^ s ( Dame et al.ll200lll to 
calculate the molecular column density based 
on the (J=l-0) emission. In the second 

method, the column density is converted to 

the H 2 c olumn density using N (ll 2 )/N(^^CO) « 
6.2 X 10^ ( Nagahama et al.]ll998ll . and in the third 
method, the relation 7 V(H 2 )/fV(C^® 0 ) ks 7 x 10 ® 
( Warin et al ][l99i) is used for the conversion from 
C^®0 column density to H 2 column density. In 


Our analysis has revealed that SNR Kes 41 is 
confined in a cavity enclosed by MCs in the west 
and the north as well as an H I cloud in the south¬ 
east (' 03.2.11) . Some shell-like features of molecular 
gas also appear to follow the periphery of the SNR. 
We can give some crude estimates on the gas den¬ 
sity of the complicated surrounding environment. 

Assuming the LOS size of the ~ —50 kms“^ 
MC is similar to its apparent size > 9' (a lower 
limit given by the FOV) and using the molecu¬ 
lar column density 7 V(H 2 ) ~ 3-5 x 10^^cm“^, we 
estimate the molecular density as n(H 2 ) < 310- 
510(i("2^ cm“®, where di 2 = (i/( 12 kpc) is the dis¬ 
tance to the MC associated with the SNR in units 
of the reference value estimated from the maser 


observation (jKoraleskv et al.lll998l) . 

Another density estimate can be made from 
the excavation of the molecular gas in the SNR 
extent , following the method used in Jiang et'HI 
( 2 OIOII . We plot the molecular column density 
distributions of the ^^CO line complex in inter¬ 
val — 70—40 kms“^ along the two diagonals of 
the FOV using X-factor method (see Figure ITUl). 
As can be clearly seen, A(H 2 ) generally increases 
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both from the southeast to the northwest and from 
the northeast to the southwest. Notably, there 
are peaks/bumps in column density just outside 
the SNR boundary, which is consistent with the 
molecular shells described in § 13.2.11 Outside the 
bumps of the both iV(H 2 ) curves, there are plat¬ 
forms at ~ 4 X lO^^cm”^; just inside the oppo¬ 
site boundaries, there are lower platforms (wide 
in the left panel and narrow in the right panel) 
both at ^ 2.9 x lO^^cm”^. The common differ¬ 
ence of column density ~ 1.1 x 10^^ cm“^ can be 
explained with the removal of molecular gas from 
the SNR region. Thus the molecular gas origi¬ 
nally in the excavated region has a number density 
n(H 2 ) ~ 140dj"2^ cm“^ assuming the mean LOS 
length of the excavated region is similar to the 
minor axis of the elongated SNR. 

If the apparent size of the southeastern H I 
cloud — that is, ~ 18' (see § 13. 2. ID — is adopted 
as the LOS size of the cloud, then the number den¬ 
sity of the neutral gas can be estimated from the 
neutral column density as n(H I)~ 36dj”2^cm“^. 
In the X-ray spectral fitting, a high hydrogen col¬ 
umn density of ~ 5.5-8.0 x lO^^cm”^ is needed 
to explain the absence of soft (< 2.0 keV) X- 
ray emission. Apparently, the H I gas to¬ 
wards the Kes 41 direction, which has a column 
density ^ 2 x 10^^ cm~ ^ I Kalberla et al. 2005; 
Dickey fc Lockmanl 1990l l H, is not sufficient to 
account for the X-ray absorption. However, the 
14 sr ~ — 50kms“^ MC, which is associated with 
the SNR and seems to dominate the CO emission 
in the direction (Figure[6]), has a hydrogen column 
density < 7-10 x 10^^ cm“^ and can significantly 
contribute to the absorbing column. Moreover, 
the harder X-ray emission in the SW than in the 
NE (§ 13.1.ID may be consistent with the heavier 
absorption by the denser molecular gas located in 
the west. 

4.2. SNR physics 

4-2.1. Comparisons with main mechanisms and 
models 

(i) Some physical parameters. The density of 
the X-ray emitting gas interior to the SNR can 
be estimated using the volume emission measure 
obtained from the spectral analysis. On the as¬ 


®See http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pi 


sumption Ue ~ 1.2nH for the densities of elec¬ 
trons and hydrogen atoms, we get ^ 0.2- 
cm“^, where / is the filling factor 
of the hot gas. We estimate the ionization age 
ti ^ Ax 10^-1.1 X 10®/^/^(ij2^ yr by using the best- 
fit ionization parameter Uet (also see Table [5]). 

In the X-ray spectral analysis, the spectral ex¬ 
traction region covers the main part of the X- 
ray emission. The mass of the hot gas in the 
region is ^ 2-4/3/2d^/2^g, which can be re¬ 
garded as a lower limit. If the hot gas pervades 
the whole interior of the remnant, we obtain an 
upper limit of the mass of the hot gas, ^ 18- 
37f/^dli^MQ, assuming a prolate ellipsoid vol¬ 
ume (12.6 X 8.4 X 8 .Adi 2 pc^) for the remnant. 

(ii) Sedov evolution. Kes 41 has been shown 
to be surrounded by molecular gas of density 
n(H 2 ) ~ 140-350cm“^ and H I gas of den¬ 
sity ^ 40cm“^. If the SNR follows 

the ISedov ( 1959ll evolution in the dense ambi¬ 
ent gas, and the hot gas temperature obtained 
from the X-ray spectral fit, kT^, is taken to be 
the emission-measure weighted mean tempera¬ 
ture, then the postshock temperature would be 
kTs = kTx/1.27 and the expansion velocity would 
be Vs = {16kTs/3p,mu)^^'^ ^ lO^krns"^, 

where the mean atomic weight p, = 0.61 is 
used. The dynamic age would he t = 2r/5vs ~ 
5 .8j'j°'3 X 10 ^di2 yr (where r ~ 13di2 pc is the mean 
radius), and the supernova explosion energy would 
be if = (25/4^)1.4nomHr^z;s ~ 2 d \2 x 10^^ erg or 
2(if 2 X lO®'^ erg, for ambient density = 40 cm“^ 
(for H I gas) and 280cm“^ (for molecular gas), 
respectively, where f = 2.026. Such estimates are 
highly unphysical, because the explosion energy 
is much higher than the canonical value 10^^ erg. 
Moreover, the molecules in the surrounding molec¬ 
ular shell (§ 13. 2. ID could not survive if the shell 
were swept up by the SNR shock, given the high 
shock velocity. Therefore, the assumption of Se¬ 
dov evolution of Kes 41 in a uniform medium is 
problematic. 

(Hi) Evolution in radiative phase. The ra- 
diatively cooled rim is one of the scenarios in¬ 
voked to explain the X-ray thermal compos¬ 
ite morphology. If Kes 41 is not considered 
to be in a pre-existing cavity, then the ra¬ 
diative pressure-driv en snowplow (PD S) stage 
begins at a radius ( Cioffi et al. Il988ll Cp^s = 
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2.9£’51^(^o/ 40 cm PC and an age 

tpDs = 3.6 X 10^£’g{^^(no/40cm“3)“'^/'^Cm®^^^yr, 
where E^i = ^^/(lO^^erg) and Cm metallicity fac¬ 
tor of order unity. If Ug = 280 cm“^ for the 
detected molecular gas, then rp^s ^ 1.3 pc and 
tpDs ~ 1.8 X 10^ yr. Compared to the mean radius 
~ 13 pc and the ionization age > 4 kyr, these es¬ 
timates indicate that the SNR may have entered 
the radiative phase at an early time. In this sce¬ 
nario, the present velocity of the blast wave is Vs = 

21 f£;si (nn/40 cm-3 ) -36/3iCm'^^^Vr/13pc) -9^/^M^^/'^^ km 
( CiofS et al. 19881 : IChen fc Slane 2001 ). which 
would be ^ 4kms“^ for an ambient molec¬ 
ular gas with Ug = 280 cm“^, where E^i ~ 

1 and Cm ~ 1 are assumed for simplicity. 

With such a shock velocity, the molecular shell 
can plausibly be explained to be the mate¬ 
rial swept up by the SNR blast wave. The 
dynamic age of the remnant is thus given by 
t = 3.3 X 10^(r/14pc)(us/413kms“^)“^[3 -|- 
nii~io/3iCm^O' ^ ^i(r/1 4 pc)-io/3i(r;,„/413knr s~^)^°/^^] yr 
( CiofS et al. 1988 : Chen &: Sland 200 ill , which 
yields t ^ 1.8 x lO^yr for a preshock atom gas 
with Ug ^ 40cm“3 or 1.0 x 10® yr for a molecu¬ 
lar gas with Ug ~ 280cm“3, respectively, where 
r = 13 pc is adopted. The age estimates here, 
however, appear to be larger than the ionization 
age, ti ~ 4-110/3/2^3/2 

(iv) Thermal conduction. Thermal conduction 
may play a role in shaping the centrally bright¬ 
ened thermal X-ray morphology of Kes 41 by low¬ 
ering the temperature and increasing the den¬ 
sity of the hot gas in the central portion. If 
the suppression of conduction by magnetic fields 
can be ignored, the conduction timescale of the 
hot gas inside the remnant is given by tcond ~ 
nekr"^/ k, where the conductivity is given hy k = 

1.8 X 10“®T,^^^/In A, with the Coulomb logarithm 
InA = 29.7-f lnn/^/^(r,,/10®K) (ISpitzerlfl^ . 

The timescale would thus be tcond 1.5-3.1 yr, 
shorter than the ionization age of the X-ray emit¬ 
ting gas. 

(v) Cloudlet evaporation. As another popular 
scenario for thermal composites, cloudlet evapo¬ 
ration is possibly taking place in the remnant in¬ 
terior, since Kes 41 is evolving in an interstellar 
environment with molecular gas which is usually 
clumpy and dense clumps can be engulfed by the 
blast wave. In this case, hot gas with a mass up 


to ^ 18-37M0 is mostly composed of evaporated 
gas. This mass is not high enough to substantially 
dilute the supernova ejecta of a few solar masses, 
and therefore the evaporation scenario may not 
be in conflict with the detected metal enrichment. 
Limited X-ray counts available in the archived ob¬ 
servations prevent an estimate of the spatial dis¬ 
tribution of such physical parameters as the tem¬ 
perature and density of the hot gas. Even the 
X-ray surface brightness shows severe asymmetry 
j. primarily due to the presence of two brightness 
peaks IS I3.1.ip . which also makes it difficult to lo¬ 
cate the explosion center. All these factors hamper 
us from presenting a quantitative comparison be¬ 
tween the observational pro perties and the evapo - 
ration effect predicted in the lWhite fc Lona ( 199ll) 
model. 

4 . 2 . 2 . Metal enrichment 

Metal enrichment in the thermal composite 
SNRs complicate the explanation of the physical 
mechanism for this class, but meanwhile it is also 
an important clue to understanding this class of 
sources. Kes 41 is a thermal composite SNR where 
the abundances of a few metal species (S and Ar) 
are elevated. This result motivates its classifica¬ 
tion as new member of the subclass o f “enhanced- 
abundance” (iLazendic fc Slanj 2006[) or “ejecta- 
dominated” (jPannuti et al.l 2014a ) thermal com¬ 
posite SNRs. This subclass has members amount¬ 
ing to a number of 25-30 (see an updated list of 
thermal composite SNRs in Table |4]) , comprising 
about 67%-81% of the 36-37 known thermal com¬ 
posites. Sulfur and silicon are the most common 
enriched species seen among these sources: they 
are overabundant in 18-19 and 14-15 such SNRs, 
respectively. The inert elments such as Ar and Ne 
are overabundant in thirteen thermal composites. 


In the case of SNR W44, IShelton et al.l l)2004al) 
suggest that central metal enrichment is a factor 
as important as thermal conduction for the expla¬ 
nation of the centrally-brightened X-ray emission 
from thermal composite SNRs. They interpret the 
enhanced metal abundance of W44 as originating 
from destroyed dust grains within the SNR as well 
as ejecta enrichment. In Kes 41, the dust destruc¬ 
tion scenario may not be favorable because argon, 
obviously enriched, is an inert element and cannot 
easily depleted into dust. In the most recent X-ray 
studies of G352.7—0.1 and W51C, metal enrich- 
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ment is explained to come from superno va ejecta 
( Pannuti et al. 2014al : ISasaki et al. 2014 k 

As pointed out by I Vink (2012). metal-rich 
plasma produces more X-ray emission than plasma 
that is not metal enhanced. In Kes 41, the S and 
Ar emissions contribute 7% and 2%, respectively, 
to the broad energy ranges 0.5-10 keV, and the 
extra abundances of these two metals contribute 
3% and 1%, respectively. 


4-2.3. Effect of pre-existing cavity 


(i) Shock reflection. Thermal composite SNRs 
are usually associated with dense material (espe¬ 
cially MCs), and therefore their progenitor stars 
are likely to have created a cavity in the dense 
environment with energetic stellar winds and ion¬ 
izing radiation. Actually, about ten thermal com¬ 
posites have been found or suggested to evolve in 
cavities (see Table 11. When the supernova blast 
wave hits the cavity wall, a reflected shock can be 
expected to be sent back inwards and reheat the 
interior gas including the ejecta to a higher tem¬ 
perature. This is the case suggested for the ther- 
mal X-ray compos ite morphology to SNR Kes 27 


(jChen et al.N2008r) . A similar explanation is also 


proposed for the thermal composite SNR 3C 397, 
which is embedded in the edge of a gi ant MC and 


confi ned by a sharp molecular wall (j.Tiang et al 
201fll l. The thermal composite SNR W44 ap¬ 


pears to be confined in a molecular cavity that 
is delineated by sharp molecular “edges” in the 
east, southeast, and southwest, as well as molecu¬ 
lar fil aments and clum ps in other periphery posi¬ 
tions ( Seta et iil]l2004h . The metal-enriched X-ray 
emitting gas in the interior of the SNR may have 
been reheated by the shock reflected from the cav¬ 
ity wall (which is at least represented by the sharp 
“edges”). 


Kes 41 has been demonstrated here to also 
be confined in a cavity in a molecular environ¬ 
ment. As shown above fS I4.2.1|) . both the Sedov 
and radiative-phase evolution in a uniform dense 
molecular/H I gas will lead to unreasonable phys¬ 
ical parameters, and hence the cavity is not blown 
by the SNR itself. A collision of the SNR shock 
with the pre-existing cavity wall thus cannot be ig¬ 
nored. The forward shock of the SNR is drastically 
decelerated upon the collision. Assuming that the 
density ratio between the cavity and the wall is 
much smaller than 1 and adopting the cavity ra¬ 


dius Tc 13 pc, the shock velocity u pon the col¬ 


lision can be essentially estimated as (jChen et al 
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( 1 ) 


which is < 20£'g(^ and SOiSj^^kms"^ for tiq > 
280 cm“^ (for molecular gas) and Uq ~ 40cm“^ 
(for H I gas), respectively. Such low velocities in¬ 
dicate that Kes 41 has entered radiative stage right 
after the blast wave struck the cavity wall. This 
effect readily explains the lack of X-ray emission 
in the SNR rim. 


The interior gas remains hot and is elevated 
to a higher temperature by the heating by a re¬ 
flected shock. This scenario can be tested in com¬ 
bination with the ob servat i on res ult and simplified 
theory presented by Sgro ( 1975ll . For a shock re¬ 
flection at a dense wall, the temperature ratio be¬ 
tween the post-reflected-shock gas and the post- 
transmitted-shock gas is given by 


Tt 



( 2 ) 


where Ar is the density contrast between the post- 
reflected-shock gas and the post-incident-shock 
gas and A is the wall-to-cavity density contrast, 
which in turn is a function of A^, specifically. 


^ __ 3Ar(4Ar 1) _ ,2, 

" {[3A,(4 - Ar)Y/^ - Vl{Ar - 1)}2 ■ 

We take n(H 2 ) 140 cm“^ as a typical den¬ 

sity of the molecular gas and adopt the velocity 
of the transmitted shock in the molecular cav¬ 
ity wall to be ^ 20kms“^, which corresponds to 
kTt ^ 2.2 eV. If the observed temperature of the 
hot gas, is adopted for T^., we can calculate 
the dependence of the wall-to-cavity density A on 
Tr in the uncertainty range 1.3-2.6 keV. FigurefTT] 
shows that A is in the range of 1500-3000. Thus 
a cavity density ^ 0.1-0.2cm“^ is implied. If 
we adopt the observed upper limit 500 cm“^ for 
n(H 2 ), almost the same cavity density is inferred. 
Such a cavity density seems to be of order similar 
to the observed hot gas density. Here, shock reflec¬ 
tion from the dense molecular wall is more efficient 
than from the H I cloud; actually, the centroids of 
both the diffuse X-ray emission (§ 13.1.11) and the 
EW maps of S and Ar (S I3.1.2I) in the southwestern 











































half can be consistent with the reheating primar¬ 
ily by the shock reflected from southwestern dense 


MC. 


(a) Progenitor’s mass. The pre-existing cavity 
of Kes 41 should be excavated by the progenitor 
star, and the progenitor’s mass can be estimated 
from the cavity size. A linear relation has been 
found between the progenitor’s mass, Mp^og, and 
the radius, Rb, of the bubble blown by the mai n- 
sequence wind in a giant MC (IChen et al.ll2013r) : 


Rb « (1.22 ± Q.05){Mprog/MQ) - (9.16 ± 1.77) pc. 

(4) 

If we adopt Rb ~ r^, the mean radius of the cavity, 
a stellar mass 18 ± 2Mq is obtained. However, 
the cavity is not purely confined by molecular gas, 
but by a giant H 1 cloud in the southeastern side, 
and there seems to be a blow-out morphology in 
the northeast (see Figure [1}. Therefore, this mass 
estimate should be regarded as a lower limit of 
the progenitor’s mass if the progenitor was not in a 
binary system, which implies that the progenitor’s 
spectral type is no later than BO. 

(iii) Ejecta in cavities. It can be seen in Ta¬ 
ble m about half out of the 36-37 known ther¬ 
mal composite SNRs are interacting with adja¬ 
cent MCs; moreover, eleven or twelve of them 
are known/thought to be evolving in cavities in 
dense gas, and most of them exhibit overabun¬ 
dance of metal species. This fact can be explained 
in the context of the ejecta in pre-existing cavi¬ 
ties. An appreciable fraction of supernova ejecta, 
especially some heavy metal species like Si, are ex¬ 
pected to expand a t a relatively low velocity (e. g., 

Mao. J. et al. 


lO^kms e. g.. lOno et al 


ei y low ■ 

DIm, 


in preparation); for the ejecta expanding in pre¬ 
existing cavities, the slow part may have not been 
substantially mixed or diluted by the interstellar 
materials by the time when the blast shock is re¬ 
flected backwards from the cavity wall. They are 
not only heated by the reverse shock in the early 
free expansion phase, but also can very possibly 
be reheated by the reflected shock. Therefore, the 
hot interior are thus observed to be metal enriched 
in X-rays. 


4 . 2 . 4 . Recombining effect? 

A number of thermal composite SNRs have 
been revealed to contain recombining plasma (see 
the rows labeled with “A” in Table H]). However, 


the over-ionization in the case of Kes 41 can be nei¬ 
ther proven nor disproven by the spectral fitting 
due to the lack of obvious radiative recombination 
continuum (RRC) features and H-like lines of the 
same metal species in the X-ray spectra with in¬ 
sufficient counts. If there is recombining plasma in 
the remnant, the temperature obtained here using 
the vnei model may be an overestimate in view 
of the “apparent” hard continuum raised by the 
RRC. Deeper X-ray observations are needed to de¬ 
termine if the X-ray-emitting plasma of Kes 41 is 
indeed overionized. 

5. Summary 

We have carried out an X-ray observational 
analysis by revisiting the archival XMM-Newton 
and Chandra data of the thermal composite 
SNR Kes 41. We have also conducted and an¬ 
alyzed an observation towards the remnant in 
^^CO, ^^CO, and C^®0 lines. The main results 
are summarized as follows. 

1. Thermal X-ray emission from Kes 41 is de¬ 
tected above 2 keV, which is characterized 
by distinct S and Ar He-like lines and a 
brightening in the SNR interior, with a cen¬ 
troid in the emission toward the southwest 
portion of the SNR. 

2. The X-ray emitting gas can be described as 
an optically thin plasma (at a temperature 
1.3-2.6 keV, with an ionization parameter 
0.1-1.2 X 10^^ cm“^ s) with a NEI model. 
The metal species S and Ar are found to 
be over-abundant, with abundances 1.2-2.7 
and 1.3-3.8 solar, respectively. They may be 
enriched by the supernova ejecta. 

3. The SNR is found to be associated with 
a giant MC at a systemic LSR velocity 
—50kms“^ and is confined in a cavity delin¬ 
eated by a northern molecular shell, western 
concave MC with also a discernable shell, 
and a southeastern H I cloud. 

4. The cavity is pre-existing prior to the super¬ 
nova explosion; it would not be physically 
realistic to conclude that it was generated 
by the SNR expansion. 

5. Kes 41 seems to have left the adiabatic 
stage of evolution and entered the radiative 
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phase with an X-ray dim rim just after the 
SNR shock encountered the dense interstel¬ 
lar medium cavity. 

6 . Thermal conduction and cloudlet evapora¬ 
tion seem to be feasible mechanisms to inter¬ 
pret the X-ray thermal composite morphol¬ 
ogy, while the scenario of gas reheating by 
shock reflected from the cavity wall is quan¬ 
titatively consistent with the observations. 

7. The birth of the SNR in a cavity in a molec¬ 
ular environment implies a mass > 18Mq for 
the progenitor if it was a single star. 
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grant BK20141044 by Jiangsu Provincial Natural 
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Table 1 

Archival Chandra observations for Kes 41 region 





Observation 

Exposure 

CCD 

Good Exposure 

Spectral 

Obs. ID 

R.A. 

Decl. 

Date 

(ksec) 

chips ^ 

(ksec) 

analysis ^ 

12513 

16''39™53?81 

-46°57'4570 

06/27/11 

20.4 

10 

20.2 

yes 

12514 

16*'39™01?38 

-46°49'4578 

06/10/11 

20.0 

11 

19.8 

no 

12516 

16''39™06!95 

-46°06'42"5 

06/11/11 

19.8 

12 

19.5 

no 

12517 

16''38™14!49 

-46°58'4270 

06/11/11 

19.8 

13 

19.5 

yes 

12519 

16*'38™19?83 

-47°15'3878 

06/13/11 

19.6 

S2 

19.3 

no 

12520 

16''37™27!34 

-47°07'3770 

06/13/11 

19.6 

S3 

19.0 

yes 


^Indicates the chips which sampled emission from Kes 41. 

^Indicates which observations featured data that were extracted for spectral analysis. 


Table 2 

Spectral fit results with 90% confidence range for the X-ray emitting gas in SNR Kes 41 


Parameter 

Value 

Net count rate (10“^ counts s) 

28.05 ± 1.60*^ 
7.86 ±0.49'" 
2.46 ±0.38® 

Vh( 10^^ cm“^) 

6.75tlll 

fcr,(keV) 

l.SOt^oir 

net(10^^s cm“^) 


/ne'RHV/dL(10^'^cm“3) 

1.46^1^2 

[S/H] 


[Ar/H] 

^•^^-1.06 

Reduced (d.o.f.) 

0.99 (54) 

Flux‘d (10“^^ erg cm“^ s“^) 

1.9 

(cm“3) 

0.3 ±0.1 

(10"^ yr) 

1 Q+8-9 
^■^-1.5 


a,&,cT]^e count rates are from PN, MOS, and 
ACIS, respectively, with Itr error range. 

‘^The unabsorbed fluxes are in the 0.5-10 keV 
band. 

®In the estimate of the densities, we assume 
a volume of a prolate ellipsoid (considering the 
NE-SW elongated morphology of the SNR) with 
size 5.8 X 4.2 x 4.2df2PC^ for the elliptical region 
of X-ray spectral extraction (see Figure |4]). 


Note. —The abundances of all elements are 
set to solar abundances except for S and Ar. 
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Table 3 

Derived Parameters of the Ambient Clouds in Velocity Interval —70—40 kms“^ 


Lines 

TpeaklK)" 

V(H2) (10^2 cm-2) 

^^(lO^Mo)^ 

T 


9.5 

3.8 

21.1d?2 

- 


4.5 

3.5 

19.5d?2 

0.64 

Ci®0 

0.95 

5.1 

N{R I) (lO^nm-'") 

28.3d?2 

0.10 

HP 

no 

5.1 

1.9df2 

- 


'‘The average peak temperature in the box region decipted in 
Figure [3 corrected to radiative temperature already. 

^di 2 = (i/(12kpc) is the distance to the MC associated with the 
SNR in units of 12 kpc. 

''We use the same LSR velocity range —70—40 kms“^ for the 
HI gas, and this range likewise includes the main part of the HI 
line profile peaked at —50kms“^. 
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Fig. 1.— Tricolor image of Kes 41 in multiwavelength. Blue: XMM-Newton X-ray image (2.0--7.2 keV) 
obtained by combining the EPIC-PN and MOS data. Green: (J=1“0) integrated map (Ksr = —70- 

—40 kms“^). Red: H I line emission from SGPS integrated map (Ksr = —55- —50 kms“^). The 
image is overlaid wi th MOST 843 MHz radio contours (at seven linear scale levels between 0.00 and 
0.79 Jy beam"^; from Whiteoak fc Green (1996)). The white cross indicates the location of 1720MHz OH 


maser ( Koraleskv et al. 1998h in this SNR.. 
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Fig. 2.— The XMM-Newton (MOS+PN) X-ray images of Kes 41 in energy bands 2.0-2.8keV, 2.8-3.5keV, 
and 3.5-7.2keV, with radio contours (the same as in Figure [1]) overlaid. 



2 5 

Energy (keV) 


Fig. 3.— X-ray spectra of the region defined in Figure|4]from both XMM-Newton and Chandra observations. 
Red points show the EPIC-PN data; the black points are from the merged EPIC-MOSl/2 data; the green 
points show the co-added spectrum made by merging all 3 Chandra observations that cover the majority 
of the SNR’s X-ray emission as described in Section 13.1.21 All the spectra have been adaptively binned to 
achieve a background-subtracted S/N (signal-to-noise ratio) of 3 per bin. They are jointly fitted using the 
same absorbed vnei model and the lines with different colors are fitted model lines matching the same color 
data points separately. 
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Fig. 4.— Raw images of the XMM-Newton observation (left) and merged Chandra observations (right), with 
all detected point sources removed. The common ellipse (in green) is used for source spectral extraction. 
In the left panel, the dashed boxes, excluding the dashed circle with the crossing line (in red), are used for 
PN and MOS1/MOS2 background subtraction. In the right panel, the boxes with the numbers denoting the 
Obs. IDs in Tabled] are used as background regions for the corresponding data. The cyan contours are from 
radio observation, with levels the same as in Figured! 
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0 0.19 0.38 0.57 0.76 0 . 9:0 0.18 0.35 0.53 0.71 0.8 

Fig. 5.— Equivalent width maps of the sulfur line (~ 2.43 keV, left panel) and the argon line 3.16 keV, 
right panel), both overlaid with the same radio contours as those in Figure [T] The energy range we use to 
extract the sulfur line emission is 2.33-2.60 keV, while the lower- and higher-energy ranges for estimating 
the continuum emission at the line energy range are 2.27-2.30 keV and 2.63-2.83 keV, respectively. For the 
argon line, the line energy range is 2.87-3.30 keV, and the lower- and higher-energy ranges are 2.63-2.83 keV 
and 3.40-3.90 keV, respectively. The color bars are in linear scale in units of keV. 
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Fig. 6.— The spectra of (J=l-0), (J=l-0) and C^®0 (J=l-0) line emission in the FOV. The 

black dotted vertical line denotes the LSR velocity of the OH maser emission. 
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Fig. 7.— (J=l-0) intensity maps integrated each 1 kms“^ in velocity range —70—40 kms“^ in linear scale, 

overlaid with the MOST 843 MHz radio continuum contours (at 7 linear scale levels between 0.07 and 0.88 Jy beam“^). 
The cross indicates the location of the OH (1720 MHz) maser spot. The colorbar of the last row of panels is different 
from those of the other rows due to the relatively low brightness. The dashed lines in the —54.5kms~^ panel depict 
the two diagonals of the FOV along which the column density distribution A^(H 2 ) of the molecular gas in velocity 
range —70—40 kms“^ are measured. The intersection of the diagonals is used as the reference point (see Figure ITOl) . 
The dashed box in the —53.5kms~^ panel decipts the region we derive the parameters of the interacting MC. 
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Fig. 8.— Intensity maps of (J=l-0) observation between —66 kms ^ and —41 kms ^ in linear scale. 
The contours and cross are the same as those in Figure [T] Note the difference in colorbars for different rows. 
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Fig. 9.— Intensity maps of C^®0 (J=l-0) observation between —55 kms ^ and —45 kms ^ in linear scale. 
The contours and cross are the same as those in Figure [71 




Fig. 10.— The column density distribution along the dashed lines illustrated in the —54.5kms“^ panel of 
Figure [71 the left panel corresponding to the northeast-southwest oriented line and the right panel to the 
southeast-northwest oriented line. The reference point is the intersection of the lines. The shaded parts 
indicates the range within the SNR’s radio shell. 
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Fig. 11.— The relation between the wall-to-cavity contrast and the temperature after the heating by the 
reflected shock. The calculation for two exemplified molecular densities is given. There are two sets of 
solutions of Equation ([3]) for each density. 
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Table 4 

List of Known thermal composite SNRs 



Enriched Metal 

Temperature 

ionization parameter 

MC 

Pre-existing 

Progenitor’s 

SNR Name 

Species 

fcTe (keV)® 

net (s cm“3)^ 

Interaction*" 

Cavity/Shells (Radius/pc) 

Mass (Mq) 

GO.0+0.0 (Sgr A East) 

S, Ar, Ca, Fe, Ni [1,2,3,4] 

fcTi ~ 1 , 
fcT2 ~5 [3,2,4] 


Y, OH 

Y (~ 5) [5] 

13-20 [1,2,3] ,> 12*^ 

G6.4-0.1 (W28) 

- 

kTi ~0.4, 

(net)i ~ 1-4 X 10^^ 

Y, OH 





kT2 ~0.8 [6,7,8] 

[6,7,8] 




R-. 

- 

~0.6 [6,9] 

~ 9 X 10^1 [6,9] 




G21.8-0.6 (Kes 69) 

- 

~ 0.8 [10,11,12] 


Y, OH 

Y C~ 13) [13] 

i) 18 [14] 

G31.9+0.0 (3C 391) 

- 

~ 0.5-0.8 [15,16,8] 

> 10^^ [15,8] 

Y, OH 

Y (~ 5) [18] 

> 12*^ 

R-. 

Mg? [17] 

~ 0.5-0.6 [17] 

~ 10^2 [17] 




G33.6+0.1 (Kes 79) 

S, Ar [19,20] 

kTi ~ 0.2, 

kT2 ~ 1.1 [20,21,19,22] 

(net)2 ~ 3.9 X 10^° 
[20,21,19,22] 

Y? 

Y (~ 8) [23,20] 

^ 14 [20] 

G34.7-0.4 (W44) 

Ne, Si, S [24] 

kTi ~ 0.3-0.8, 

kT2 ~ 0.9-4 [24,8,28,25,27] 

> 3 X 10-*-^ [24,8,25] 

6 X 10^^ [29] 

Y, OH 

Y (~ 30) [26] 

8-15? [27] 

R-. 

Mg, Si, S, Ar, Ca [29] 

kTi ~ 0.5, 
kT2 ~ 1.0 [29] 





G38.7-1.4 

O, Ne [30] 

~ 0.6 [30] 

~ 7 X 10^^ [30] 




G41.1-0.3 (3C 397) 

O, S, Ca, Fe [31,32,33,34,35] 

kTi ~ 0.21-0.25, 

(raet)l ~ 3 X 10^^, 

Y 

Y (~ 6) [36] 

> 12 [14] 



kT2 ~ 1.5-3.5 [31,35,34] 

(net)2 ~ 2 X 10^^ [31,35,34] 




G43.3-0.2 (W49B) 

Si, S, Ar, Ca, Fe, Ni [8,37,38,39,40] 

kTi ~ 0.7-1.05, 
kT2 ~ 1.75-3.3 
[39,41,42,38,8,40] 


Y? 

Y? (~ 7) [23] 

> 13 [23] 

R: 

Si, S, Ar, Ca, Fe, Ni [43,44,45] 

kTi ~ 0.13-0.3, 

kT2 ~ 1.12-1.91 [43,44,45] 




~ 25 [37,39,41,40] 

G49.2-0.7 (W51C) 

Ne [46], Mg [46,47], S [48] 

0.56-0.74 [46,47,49,48] 

0.8-11x10^^ 

Y, OH 


> 20 [46] 

G53.6-2.2 (3C 400.2) 

Fe? [50] 

~ 0.8 [50] 

~ 10^^ [50] 




G65.3+5.7 

- 

~ 0.25 [51] 





G82.2 + 5.3 (W63) 

Mg, Si, Fe [52] 

kTi ~ 0.2, 
kT2 ~ 0.6 [52,53] 





G85.4+0.7 

O? [54] 

1.0 [54] 

8X10^*^ [54] 




G85.9-0.6 

O, Fe [54] 

1.6 [54] 

5x10-^^ [54] 




G89.0+4.7 (HB21) 

Si, S [55,56] 

0.62-0.68 [56,55] 

> 4 X 10^^ [55,56] 

Y 



G93.3+6.9 (DA 530) 

Si [57] 

0.3-0.6 [57,58] 

> 4 X 10^^ [57] 



~ 10 [57] 

G93.7-0.2 (CTB 104A)? [53] 

- 






G116.9+0.2 (CTB 1) 

Mg [55], O, Ne [56] 

kTi ~ 0.2-0.3, 
kT2 ~ 3 [56], 
or kT 2 ~ 0.8 [55] 

> 1 X 10^^ [56,55] 



> 13 [56] 

G132.7+1.3 (HB 3) 

Mg [55], or Mg, Ne, O [55] 

~ 0.3 [55, 53] 


Y? 

Y (~ 70) [59] 


G156.2 + 5.7 

Si, S [60,61,62] 

kTi ~ 0.45, 

kT2 ~ 0.6 [62,60,63,61] 

~ 1.5 X 10-^-^ [60,61,62] 



< 15 [61] 

G160.9+2.6 (HB 9) 

- 

~ 0.8 [64] ^ 


7 



G166.0+4.3 (VRO 42.05.01) 

S [65] 

~ 0.7 [66,67] 


7 



G189.1 + 3.0 (IC443) 

Mg , S [68,69,65], Si [68,69], Ne [65] 

kTi ~ 0.3-0.7, 
kT2 ~ 1.0-1.8 
[68,69,65,70,71,53] 

(net)i > 10^"^ [68,69] 

~ 9.8 X 10^1 [73] 

Y, OH 

Y (~ 11) [75] 

> 15 [75] 

R-. 

Si, S, Ar [72], Ca, Fe, Ni [73] 

~ 0.65 [72,73,74] 




G272.2-3.2 

O [76], Ne [77], Si, S, Fe [76,77,78], Ca, Ni [78] 

~ 0.7-1.5 [77,76,78,79,80] 

2-10x10-^'^ [77,76,78,80] 



(SN la) [78,77,76] 

G290.1-0.8 (MSH 11-61A) 

Si, S [81,82], Mg [81] 

~ 0.6-0.9 [82,81] 

> 0.1 X 10^^ [82] 

7 


25-30 [82] 

G304.6+0.1 (Kes 17) 

Mg? [83,84,85] 

~ 0.79-1.10 [84,86,85,83] 

> 3.7 X 10^^ [84,86,85,83] 

Y 



G311.5-0.3 

- 

~ 0.98 [84] 


7 



G327.4+0.4 (Kes 27) 

S, Ca [87], Si [8] 

~ 0.5-1.2 [87,53,88,8,89] 

> 1 X 10-^-^ [87,8] 


Y [90] 


G337.8-0.1 (Kes 41) 

S, Ar [91] 

~1.9 [92,91] 

~ 2 X 10-^-^ [92,91] 

Y, OH 

Y (~ 13) [91] 

> 18*^ 

G344.7-0.1 

Al, Si, S, Ar, Ca, Fe [93,94,95,96] 

~0.8-1.8 [96,95,94,93] 

1-4x10^^ [96,95,94,93] 

7 


(SN la) [93] 

G346.6-0.2 

Ca [97] 

~ 1.2 [97,98,84] 

~ 2.9 X 10^^ [97,84] 

Y, OH 



R; 

- 

~0.3 [99] 

~ 4.8 X 10^1 [99] 




G348.5 + 0.1 (CTB 37A) 

- 

0.55-0.83 [84,98,100,101,102] 

> 3 X lO-^'^ [84,102] 

Y, OH 

? [103] 


R; 

Si? [102] 

~0.5 [102] 

~ 1.3 X 10^2 [102] 




G352.7-0.1 

Si, S [104,105,19], Ar? [19,104] 

0.8-2.1 [104,105,4] 

1-2x10^'^ [104,105,19] 




G355.6-0.0 

Si, S, Ar, Ca |1061 

~0.6 [106,98] 





G357.7-0.1 (Tornado) 

- 

~0.6 [11,107] 


Y,OH 



G359.1-0.5 

Si, S [108,109] 

kTi ~ 1.1, 

kT2 ~ 2.0 [109,111,108,112] 


Y, OH 

Y 28) [llOJ 


R; 

Si, Mg, S [113] 

~ 0.3 [113] 














































Note. — “Y” in the table means “yes”; “?” means that the property is uncertain or controversial; “i?” denote that the parameters are derived using the recombining plasma model for the same SNR. 

®'The different tmperatures of different regions are shown as a range, and fcT 2 ” are temperatures of cold and hot components respectively. 

^The ionization parameter of underionized state and overionized state respectively. 

‘^This column is adopted from Jiang et al.’s (2010) SNR-MC association table. “OH” means that the 1720 MHz OH maser is detected. 

^Estimated using the method described in lOhen et al ]Gni3)- This method is not applied to SNRs W44, HB 3, and G359.1 — 0.5 which are also in molecular cavities, considering the cavity sizes are larger than the 
scope of application. 
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